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Abstract. We employ the one-dimensional Fourier representation (1DFR) to analyze 
the 3-year WMAP de-biased internal linear combination (DILC) map and its possible 
contamination by galactic foregrounds. The 1DFR is a representation of the spherical 
harmonic coefficients for each ^-mode using an inverse Fourier transform into one- 
dimensional curves. Based on the a priori assumption that the CMB signal should 
be statistically independent of, and consequently have no significant correlation with, 
any foregrounds, we cross-correlate the 1DFR curves of 2 < £ < 10 modes, which 
are claimed by the WMAP team to be free of contamination and suitable for whole 
sky analysis. We find that 8 out of the 9 modes are negatively cross correlated with 
the foreground maps, an event which has a probability of only 9/512 ~ 0.0176 for 
uncorrelated signals. Furthermore, the local extrema of the 1DFR curves between the 
DILC and those of the foregrounds for I = 2 and 6 are correlated with significance 
level below 0.04. We also discuss the minimum variance optimization method and 
use the properties of the measured cross-correlation to estimate the possible level of 
contamination present in the DILC map. 
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1. Introduction 

The cosmic microwave background (CMB) radiation contains a wealth of information 
about our Universe. Not only does the angular power spectrum of the CMB 
temperature fluctuations allow us to determine cosmological parameters that shape our 
understanding of the Universe, but the issue of primordial Gaussianity and statistical 
isotropy of the CMB also taps the most fundamental principle in cosmology, and 
any robust detection of a deviation would have far-reaching and indeed revolutionary 
implications. 

Interpreting the term "Gaussianity" in the most relevant way for cosmology, i.e. as 
meaning that the fluctuations form a statistically homogeneous and isotropic Gaussian 
random field [3, 7], the issue of non-Gaussianity in the CMB data was first raised by [29] 
in the COBE data [56] . After the release of the 1-year Wilkinson Microwave Anisotropy 
Probe (WMAP ) data [4, 5, 41, 37], departures from Gaussianity have been detected with 
various methods [12, 32, 14, 53, 28, 62, 9, 35, 47, 43, 49, 45, 61, 22]. Although Gaussianity 
(as we define it) requires statistical isotropy, a special focus and results on the breaking 
of large-scale statistical isotropy have been reported in [60, 19, 27, 15, 55, 42, 6, 1]. 
Among them, the most notable anomaly involves the alignment of the quadrupole and 
octupole alignment in the general direction of Virgo. In the 3-year WMAP data [38, 57] 
the previously detected departures from Gaussianity and statistical isotropy still persist 
[11, 17, 46, 16]. The question is whether they are attributed to foreground residuals [2], 
systematic errors [8], local structures [39] or even new physics [40, 33]. 

However, it is important to eliminate the least exciting possibilities before getting 
carried away by the most exotic ones. In particular, one has to be cautious about 
foreground residuals as they are surely present in the derived CMB map. The elimination 
of foregrounds and the extraction of CMB signal is based on the concept of minimum 
variance optimization with the a priori assumption that CMB and the foregrounds 
are statistically independent: linearly combining all available maps and minimizing the 
variance of the combined map to reduce the foreground residuals as much as possible. 
The 3-year de-biased internal linear combination (DILC) map is based on such a concept 
by the WMAP science team. This approach can certainly help to achieve the optimal 
estimate of the angular power spectrum for the CMB. It cannot, however, guarantee 
the optimal morphology (spatial distribution). Since most of the anomalies, with the 
exception of the reported low quadrupole power [24, 19, 25], are related to the pattern 
of fluctuations rather than simply their amplitude, one needs to check the derived CMB 
map very carefully before any scientific conclusion is reached. In this light, several 
authors have investigated possible foreground contaminations [48, 20, 21, 52] and the 
effect of different subtraction methods [18, 54]. 

The a priori assumption that the CMB and the foregrounds are statistically 
independent and thus have no significant cross-correlation is not only the backbone of 
the minimum variance optimization method, it is also the most basic principle (preceding 
Gaussianity of the CMB), furnishing the most fundamental statistical test for foreground 
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cleaning. Note that statistical independence involves an ensemble of universes, there 
will, therefore, inevitably be accidental cross-correlations caused by chance alignments, 
particularly on large angular scales due to the cosmic variance effect [63]. 

Efforts on the investigation on large-scale foreground contamination by direct cross- 
correlation were made in the harmonic domain. Based on the close connection between 
the phases of the spherical harmonic modes, and morphology [fO], [50, 5f, 13] examine 
foreground residuals by cross-correlation of phases between the WMAP internal linear 
combination maps and the derived foreground maps. 

In this paper, we use a new representation of the spherical harmonic modes to 
analyze the cross correlation between the DILC and foreground maps. Instead of using 
the the two-dimensional spheres, we use an inverse Fourier transform on the spherical 
harmonic modes for each I producing one-dimensional curves for each harmonic scale. 
We therefore call such representation one-dimensional Fourier representation (1DFR). 
Chiang and Naselsky [13] first devised the 1DFR to illustrate connection between phase 
coupling and morphology, particularly relating to local extrema higher than 3a. Chiang, 
Naselsky and Coles [11] use the same representation to demonstrate the anomalies in 
the distribution of global extrema of the 1DFR curves of the DILC map for the i < 10 
modes. 

This paper is arranged as follows. In Section 2 we summarize the Gaussian random 
hypothesis of the CMB. In Section 3 we introduce the 1DFR and its properties and 
advantages. We then use simple cross-correlation and extremum correlation to examine 
foreground contamination in the DILC map in Section 4. We discuss the minimum 
variance optimization method and use cross-correlation coefficients to estimate the level 
of foreground contamination in Section 5. The conclusions are presented in Section 6. 

2. Gaussian Random hypothesis of the CMB 

The statistical characterization of CMB temperature fluctuations (against the CMB 
isotropic temperature T = 2.725 K) on a sphere can be expressed as a sum over spherical 
harmonics: 



e=o m =-e 

where the Yg m (9, tp) are spherical harmonic functions, defined in terms of the Legendre 
polynomials Pe m using 



and the ae m are complex coefficients which can be expressed with ae m = \a-e m \ exp(i<3>£ m ) 
and <&i m are the phases. We use the Condon-Shortly definition for the spherical 
harmonic decomposition. Isotropic Gaussian random CMB temperature fluctuations 
on a sphere, of the type that result from the simplest versions of the inflation paradigm, 
possess spherical harmonic coefficients ag m whose real and imaginary parts are mutually 



(2.1) 




(2.2) 
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independent and both Gaussian, or equivalently, the amplitudes \a,£ m \ are Rayleigh 
distributed with random phases [3, 7]. The statistical properties are then completely 
specified by the second-order statistics, the angular power spectrum Ce, 

( a em a *e'm') = C i 5 U' 5 mm'- ( 2 - 3 ) 

Since T is always real, the complex vectors of the a^ m on the Argand plane for m < 
are mirror images of m > with respect to x axis for even m, and with respect to y 
axis for odd m. The statistics for real signal on a sphere are therefore registered only 
in the spherical harmonic coefficients ae m for m > 0. 



3. One dimensional Fourier representation of the spherical harmonic 
coefficients 

Recent studies of the large angular scale properties for the CMB temperature anisotropy 
[60, 19, 55, 42] are based on the composite map constructed for each £ by summing all 
the m modes pertaining to that £: 

£ 

Ti(0,(f)= a em Y em (6,(p). (3.1) 
m=-e 

Alternatively, one can represent a em in each £ by an inverse Fourier transform, as the 
ap m is now a function of a single variable m. Such a 1DFR for ae m from the spherical 
harmonic decomposition of the sky is written as 

i i 
Ti{<p) = a£ m exp(?my?) = a eo + 2 ^ cos(nup + $^ m ), (3.2) 

m=—£ m=l 

where for the negative m we use the complex conjugate (a£- m = a} m ) to make T^ip) 
real. As we have mentioned for real T on a sphere the statistics are registered only 
in m > modes, the T^tp) curves assembled in this way contain the same amount of 
information as the spherical Ti(0, ip). The variance for each curve is then 

°t = T- \T^)-T t ] 2 d V = 2 M 2 - (3.3) 

For a Gaussian Random Field (GRF), the \a tm \ have a Rayleigh distribution, so a\ has 
a chi-square distribution. 

3.1. Comparison between 1DFR and integration on 9 of the composite maps 

Since the 1DFR results in AT as a single function of tp, such representation is therefore 
similar to the integration on 6 used in the composite maps: 

rT t {e,ip)M. (3.4) 
Jo 
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There are, however, some subtleties. Equation (3.1) can be written as follows 

l2£ + 1 



/ 2i + 1 f 

T e (9, <p) = W | a OT ^o(cos 0) 



+ 2E 



m=l 



(-1)' 



\ TIT, V\ l a **l cos(mip + $* m )P< m (cos0) 

\ (/ + my. 



(3.5) 



As / Pe rn (cos9)d9 = for odd £ + m, the curves from the integration over from the 
composite maps have contributions from only the odd m part for odd £ and the even m 
part for even I. Therefore for odd £, the (— l) m term inverts the resulting temperature 
T as seen from the 1DFR. (As we use the Condon-Shortley definition for the spherical 
harmonic functions, / P £m {cos9)d9 is always positive for m > and for m = 0). 
Accordingly, in Figure 1, for the odd £ curves the amplitudes are multiplied by —1 to 
facilitate direct comparison. The ip coordinate of the 1DFR is also plotted backwards in 
order to match the Galactic longitude coordinate / (not to be confused with the spherical 
harmonic mode £; see Figure 2 (top panel). 

Moreover, the integration over 9 results in more symmetric curves than the simple 
1DFR, where equal footing is given to each ai m in the summation. As shown in Figure 1, 
the curves (red) obtained by integration of 9 from the composite maps show repetitions 
for up at [0,180°] and [-180°, 0]. 

3.2. Characteristics of the 1DFR 

The 1DFR assembles spherical harmonic coefficients into one-dimensional curves, so the 
information on 9 direction is lost, compared with standard two-dimensional spherical 
maps. Nevertheless, the morphology and statistics registered in the complex a £m 
sequence should still manifest themselves in this representation. For a Gaussian 
random field on a sphere, the real and imaginary parts of the spherical harmonic 
coefficients a im in each £ are mutually independent and both Gaussian distributed 
with zero mean and variance Cg/2, where Cn is its angular power spectrum. Thus 
if the a£ m are a result of Gaussian process, the 1DFR T e ((p) curves shall possess all the 
usual Gaussian random properties such as two-point correlation, peak statistics, and 
Minkowski functionals. . . etc.. 

One of the 1DFR characteristics is that the a«) modes contribute to the ID curves 
only a constant shift without altering the morphology, whereas in spherical harmonic 
composition it produces modulation in 9 direction, a£ Pio(cos 9) . This is useful when 
one is to compare two sets of ae m with standard cross-correlation, particularly when 
one of them is from a foreground map. Standard analysis on large angular scale (low 
multipole) anomalies is performed on a composite map by Equation (3.1): a full-sky 
map synthesized from the a £m from — m to m for each £ (see Figure 2 top panel as an 
example). For foreground maps defined in Galactic coordinate system, the emission near 
the Galactic plane dominates the signal, which is then spherical-harmonic transformed 
into high amplitudes of a^ for even £ (Figure 2 bottom), which produce a prominent belt 
in the composite maps, making them rather difficult to analyze under such situation. 
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Figure 1. Comparison of the 1DFR of the DILC map (blue curves) and the T 
distribution curves from integration of Ti{9, <p) on 9 direction (red curves). The x axis 
is ip and is plotted reversely to follow the conventional Galactic longitude coordinate I 
and the y axis is in unit of thermodynamic temperature ([iK). Note that the amplitude 
of the 1DFR curves for odd I is multiplied by —1 in order for comparison with those 
from integration on 8. 

4. 1DFR and Foreground contamination in the DILC map 

In this section we investigate foreground contamination in the WMAP 3-year DILC 
map. It is based on the a priori assumption that the CMB (at the background) should 
have no "knowledge" about what the foregrounds look like. So if the derived DILC map 
more or less reflects the morphology of true CMB, the DILC and foregrounds should 
have little or no resemblance to it; we can also see how clean is the derived DILC by 
directly comparing the morphology. 

Firstly, in the top mosaic of Figure 3 we plot the 1DFR curves for the DILC (blue 
curves) and those of the WMAP derived foreground maps at Q (orange), V (green) and 
W (red) channels. The foreground maps are the sum of the synchrotron, free-free and 
dust templates obtained via Maximum Entropy Method [38]. As one can see, the most 
striking feature is the significant anti-correlation between the DILC and foregrounds for 
the quadrupole (£ = 2). 

In the bottom mosaic of Figure 3 we plot, on a unit circle for each 1DFR curve, the 
positions of local extrema: all <p for dT^/dip = and d 2 T^/dip 2 ^ 0. The angle between 
each sign and the positive x axis is <p for each local extremum position. The blue filled 
signs are from the DILC and the open signs the foreground map at Q (orange), V (green) 
and W (red) channels. We use circles and diamonds to denote peaks (d 2 Ti/d<p 2 < 0) 
and troughs [d 2 Ti/dip 2 > 0), respectively. The size of the signs is ~ 5° so that one can 
see there are 12 out of totally 16 local extrema for I = 10 DILC 1DFR curves with W 
channel foreground map 1DFR local extrema located within 5°. 
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Figure 2. The composite map of t = 4 from DILC map (top) and that of I = 4 
from the W channel foreground map (bottom). One can see the dominant stripes of 
the foreground composite map caused by aeoPio{cos9), which makes it difficult to be 
compared with the DILC composite map. 



In what follows, we employ cross correlation and correlation of extrema of the 
1DFR curves to illustrate the foreground contamination in the DILC map. The cross 
correlation is a measure of "trend": two curves will yield strong cross correlation if 
both go up or down in tandem at the same interval of tp. However, any correlation will 
be cancelled out if they match on one half of ip, but on the other they are opposite. 
Although cross and extremum correlation are not totally independent, it is therefore 
intuitively helpful to examine both properties. 



4-1. Cross correlations of the DILC map with the WMAP derived foreground maps 



We use cross-correlation to quantify the foreground contamination in the DILC map: 



XI 



(27r)- 1 r 27r T&ip) - T\ Tl{p)-Tl dp 



<J,<r 



3 



(4-1) 



where Tg((p) indicates the 1DFR curve of multipole number i at p of % map, o~\ is the 
standard deviation of the curve. The X coefficients range from —1 to 1. In Figure 4 we 
show the X coefficients between 1DFR curves of the DILC and those of the foregrounds. 
First of all, if the DILC map is uncorrelated with the foregrounds, the probabilities for 
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Figure 3. Top mosaic: the 1DFR curves for DILC (blue curves) and those of the 
foreground (sum of the synchrotron, free-free and dust templates) maps at Q channel 
(orange curves), V channel (green curves) and W channel (red curves). The x axis 
is ip and is plotted reversely to follow the conventional Galactic longitude coordinate 
I and the y axis is in unit of thermodynamic temperature (/iK). Notice strong anti 
cross correlation for the quadrupole. Bottom mosaic: the position of the local extrema 
plotted on unit circles. The angle between each point and the positive x axis is the 
position <p for each extrema of the 1DFR curves from the DILC map (blue filled sign) 
and those from foreground map at Q channel (orange open sign), V (green open sign) 
and W channel (red open sign). The circle and diamond signs represent peaks (local 
maxima) and troughs (local minima), respectively. The signs subtend about 5° so 
that, for example, one can see for i — 10 there are 12 out of totally 16 local extrema 
of DILC 1DFR curves located within 5° of those of W channel foreground curve. 
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Figure 4. Cross correlation of the 1DFR curves from the DILC map with those 
from WMAP foreground maps (top panel) at K channel (green x ) , Ka channel (red 
+ signs), Q (blue A signs), V (black (}), W (black line) and cross correlation of the 
1DFR curves between the foregrounds (bottom 4 panels). Note that 8 out of 9 modes 
are negatively correlated, the probability for which to happen is Cf(2)~ 9 ~ 0.0176. 
The yellow, light blue and dark blue areas denote 1, 2 and 3 a respectively (68.27%, 
95.45% and 99.73%) of 10 5 Monte Carlo simulations. 
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positive and negative X should be equally 1/2. One can see that eight out of the nine 
modes (except £ = 10) X < 0, an outcome for which the probability is Cf(2)~ 9 ~ 0.0176. 

Secondly, the cross correlation in quadrupole is below —0.8 (at around 90% CL). In 
our previous publication [11], we use cross correlation of phases between the DILC and 
foreground maps to examine foreground contamination. The most prominent correlation 
appears in octupole, which renders the significance level as low as 0.06. Although we do 
not claim that the quadrupole-octupole alignment is directly caused by the foregrounds, 
we can at least claim that since there is significant foreground contamination, such 
alignment is not cosmological. In the lower 4 panels of Figure 4 we also show the cross 
correlation between the foregrounds, which are all above 0.96. 

4-2. Correlation of extrema of the 1DFR curves 

In our previous publication [11], we used a 1DFR on the 9 modes 2 < £ < 10 of the 
DILC map and examined the positions of the global extrema of the 9 curves. If the 
signal were Gaussian, the global extrema should distribute randomly on the (p axis. We 
found that the extrema are concentrated around ip = 180° and avoid the tp = (Galactic 
centre) with significance level below 0.5%, which is a strong indication that the global 
extrema of the DILC map are non-randomly distributed due to some influence related 
to a Galactic-coordinate frame. Here we further examine the local extrema of the 1DFR 
curves of the DILC map and their correlation with those of the foregrounds. 

In order to extract information about correlation of extrema between DILC and 
the foregrounds, we use the following method: for each foreground extremum, we search 
for its nearest DILC extremum and denote their distance by |A</^(i)|, where % denotes 
the ith extremum of the foreground curve. The collection of the separation angles for 
each £ is to be compared with the distribution of 10 5 Monte Carlo simulation on the 
GRF extrema against the foreground ones. In the top mosaic of Figure 5 we show the 
cumulative probability distribution of the separation angles between the DILC and W 
channel foreground extrema (red curves) and that of separation angles of 10 5 realizations 
of simulated GRF and W channel foreground extrema (black) . The bottom panel is the 
significance level from Kolmogorov-Smirnov (K-S) test for the distribution of extrema 
of DILC with foreground maps at K channel (green x signs), Ka channel (red + signs), 
Q (blue A signs), V (black <0> signs), W (black line). One can see that for £ = 2 and 6, 
the significance levels are all below 0.1 for all channels, and specifically, the significance 
levels are below 0.03 for the V channel, which is also the channel with the highest weight 
assigned by WMAP in the internal linear combination of multi-frequency cleaning to 
extract the CMB signal. 

4-3. Cross correlation with foreground templates independent from WMAP data 

To further certify foreground contamination in the DILC map, we also test various other 
full-sky foreground templates that are measured, collected or derived independently from 
the WMAP data. The existing foreground templates are the 408 MHz radio continuum 



The 1DFR and Large Angular Scale Foreground Contamination in the WMAP3 Datall 




20 40 60 BO 100 
[Alfrl (degrees) 




10 20 30 40 

lAi^l (degrees) 




5 10 15 20 25 
lA-^1 (degrees) 




20 40 60 

I AV I (degrees) 




10 20 30 

lA*^] (degrees) 




5 10 15 20 
lA^I (degrees) 




10 20 30 40 50 
lAV ' (degrees) 




5 10 15 20 25 30 
lA^I (degrees) 




5 10 15 20 

lA^I (degrees) 



1 .00 r 



a 



0.10 



0.01 




2 4 6 



1 



Figure 5. Top mosaic: the cumulative probability distribution of the separation angles 
(the angle of the nearest extremum) between the DILC and the W channel foreground 
cxtrcma is plotted in red curves and the cumulative probability distribution of the 
separation angles from 10 5 Monte Carlo simulation is in black. The bottom panel 
is the significance level from K-S test for the distribution of extrema of DILC with 
foreground maps at K channel (green x signs), Ka channel (red + signs), Q (blue A 
signs), V (black (} signs), W (black line). For I — 2 and 6, the significance levels are 
below 0.04 for the V channel, which is also the most weighted channel by WMAP in 
internal linear combination of multi-frequency cleaning to extract the CMB signal. 
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Figure 6. Cross correlation of the 1DFR curves from the DILC map with those from 
full-sky foreground templates that are observed, collected or derived independently 
from the WMAP data. These independent templates are the 408 MHz Continuum 
Survey (green x) for synchrotron emission, Ha map for free-free emission (red + 
signs), and the dust template extrapolated from COBE/DIRBE data (blue A signs). 
For the syncrotron and dust templates, 8 of the 9 modes are negatively correlated with 
the DILC, the probability is C?(2)- 9 ~ 0.0176. 



survey [36] tracing the synchrotron pattern, the Ha template for thermal bremsstrahlung 
(free- free emission) [30], and the dust template extrapolated from COBE/DIRBE data 
to microwave frequencies [31]. We repeat the cross correlation of the 1DFR curves 
in Section 4.1 and show the cross correlation coefficient X in Fig.6. For both the 
synchrotron and dust templates, 8 of the 9 modes are negatively correlated with the 
DILC, the probability of this happening in the absence of real correlation only being 
Cf(2)~ 9 ~ 0.0176. We find no significant evidence of cross-correlation with the Ha 
template, but this is not unexpected since this makes only a small contribution to 
the expected foreground at relevant frequencies. One should note that, although 
the 408 MHz map and the dust template are templates relating to single foreground 
components, they can be considered as total foreground templates as they represent the 
dominant component at the corresponding frequencies. Their correlation with the DILC 
strengthens the claims we made in the previous section. 

5. Cross correlations with the foregrounds as an estimate on the 
foreground contamination 

5.1. Minimum variance optimization and foreground residuals 

In the last section, we have shown through strong cross correlation and peak correlation 
that the DILC signal is probably contaminated with foregrounds. We now examine the 
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concept of multi-frequency cleaning method and how foreground residuals propagate 
into the final map. Multi-frequency cleaning has been the workhorse for retrieving the 
CMB signal from polluted data. It is based on the concept that the CMB signal exists, 
among different frequency maps, as a frequency independent component [5, 38] (see also 
[26]): 

T i = T cab + F i , (5.1) 

where T % represents the total measured signal at frequency band i, T cmh the frequency 
independent (i.e. the CMB) signal, and F l the foreground at frequency band % (here we 
assume noise is not important). This frequency-independent component can be flushed 
out with an internal linear combination on M frequency bands of maps with weighting 
coefficients Y^Li w % — 1 ; 

M M 

T ilc = = T cmb + WiF\ (5.2) 

i=l i=l 

where T llc represents the ILC map. The a priori assumption that the frequency 
independent component should be statistically independent with others ensures the 
variance of the ILC map is the sum of the variances of the two parts: 

M 

var[T ilc ] = var[T cmb ] + var[]T Wi F% (5.3) 

i=i 

Therefore if one is to minimize the variance of the ILC map by dvax[T llc ]/dwi = 0, it is 
equivalent to minimizing the variance of linear combination of the foregrounds, 

<9var[T ilc ] _ gvar[gg^Fj 
dwi dwi 

thus reducing most the residuals to produce a map closest to the CMB. The 
WMAP DILC map is thus produced by employing the minimum variance optimization 
in the pixel domain in 12 separate regions, where Region marks the largest region 
for |6| > 15°, and Region 1 — 11 for those around the Galactic plane. In each region 
the 5 frequency maps are linearly combined and a set of weighting coefficients w\ are 
obtained in such a way that the resultant variance is minimum. 

It is worth noting that Equation (5.3) requires statistical independence that is 
defined between the foregrounds and an ensemble of universes, but the fact that we have 
only one universe shall introduce some error in the minimization process. Furthermore, 
since the minimum variance optimization takes on the overall variance of the map, 
instead of the individual multipoles as the method by [59, 60, 18], there is no guarantee 
that the contamination from the foregrounds in each multipole I will be minimum (let 
alone zero). A spherical harmonic decomposition of Equation (5.2) gives 

i 

where f\ m represents the spherical harmonic coefficients for the foreground map F % . For 
higher multipoles, the power spectrum Ci of the DILC should be close to that of CMB 
as there are more m modes participating in the summation (2£ + 1) -1 J2 m \ a e m c \ 2 , as l° n g 
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as the foreground residual w ifl m ls more or l ess non-correlated between each m of the 
same £ (like white noise). On the other hand, for lower multipoles, the discrepancy 
between CMB and the DILC can be prominent when the foreground residuals are 
correlated. 

From the significant cross and extremum correlation as shown in Figure 4 and 
Figure 5, it is indeed likely that the foreground residual is still a non-negligible part for 
low multipoles in the DILC map, unless the true CMB signal happens to 'resemble' the 
foregrounds, an eventuality which has a small probability. 



5.2. Estimate of the foreground contamination 

To make an estimate of the foreground contamination present in the DILC map, we can 
first of all write 

fim = vfftm, (5-6) 

where o~j is the r.m.s. of the fluctuation of f\ m and fe m is the unitary foreground signal 
(related to morphology) with its r.m.s. of fluctuation a{ — 1. Such a degeneracy can 
be seen in Figure 4 that most cross-correlations between foregrounds are above 0.99 
(except with the W channel ~ 0.97). 

We now calculate the cross correlation coefficient of the 1DFR curves between the 
true CMB and the unitary foreground signal: 



2tt Jo o-f^al * 

2ELil« b H/^|cos(^- b -$L) 



_cmb 



(5.7) 



Although the CMB and the foregrounds should be statistically independent, i.e. 
{X c /) = 0, the fact that we only have one realization of the CMB fluctuation makes 
this parameter crucial, especially for low multipoles. 

To relate the to the measured , we first show, in the top panel of Figure 7, 
the cross-correlation of the foreground residuals with the foregrounds. The DILC map 

( Ft) 

is processed with 12 sets of weighting coefficients w\ , where i = 1 — 5 represent K, Ka, 
Q, V and W bands, respectively and (R) indicates Regions — 11, so for the foreground 
residual we use a linear combination of the WMAP foregrounds on 12 separated regions 
to obtain the whole sky map: 

my?) = EE^%), ( 5 - 8 ) 

R=0 i=l 

from which we produce the 1DFR curves Re- One can see that for £ = 2, 3 and 4, the 
foreground residual has significant anti-correlation with the foregrounds. We caution 
that the WMAP foregrounds are retrieved via the Maximum Entropy Method, not by 
direct subtraction of the DILC signal off the frequency maps, so the DILC is not a 
simple sum of the true CMB and the foreground residual as Equation (5.3). For the 
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modes in which the residual has significant correlation, we express the residuals in terms 
of the unitary foreground signals Q/g m : 

aZ c - aZ h + eefem, (5.9) 

where q is the contamination parameter at multipole number £. For other multipoles, 
this approach means we are estimating the most possible contamination. 

The cross correlation between the 1DFR curves of the DILC and the foreground 
maps can be written as 

vdf i 1 [ 2 * E m ^o a tm exp(irmp) £ m < y f , exp(imV) , 
X i =WZ I ^ n d( P 



2tt Jo o-f c af 

ee + 2 Ej = i I \hm\ cos($^ b - $L) (5 1Q) 

U\ + K mb ) 2 + 4q EL=i fe b | |/, m | cos($^ b - ${ m ) ' 

We can now denote = as it is independent of the individual foreground 

channels. The measured X t , as shown in the top panel of Figure 4, are indeed 
degenerate due to the resemblance in morphology between the foregrounds. With 
Equation (5.7) we can solve Equation (5.10) for ef. 



c — ,cmb [ y-cf . ydf 
e i - °i I ~ A £ + A <? 



(X, 



c f\2 



\ 1 - (X, 



V\2 



(5.11) 



Equation (5.11) describes the contamination parameter in relation to X'f and Xf*. 
In the middle panel of Figure 7 we plot the foreground contamination parameter ee in 
unit of <7£ mb against the cross correlation coefficients Xf . The thick line is for zero 
correlation between CMB and the foregrounds: X^ = 0, which is the case for high 
£ with more m modes in the summation of Equation (5.7). For low £, however, the 
X c / plays a major role. In the bottom panel of Figure 7 we use the mean value of 
the Xf r from DILC and WMAP foreground maps in Figure 4 as input, and run 10 
Monte Carlo simulations for the X c / . The foreground contamination is illustrated by 
the ratio e^/af^ at the values of 1, 2 and 3 a thresholds (yellow, light blue and dark 
blue areas, respectively). Note that the foreground contamination shown in this way 
is independent of the underlying CMB power spectrum, and we assume the foreground 
residuals in the DILC map has the same morphology as the foregrounds themselves, so 
we are estimating an upper limit on the contamination ratio. 



6. Conclusion and Discussion 



In this paper we have proceeded from the basic assumption that the CMB should 
have no or little correlation (within cosmic variance limit) with the foregrounds, to 
examine foreground contamination in the WMAP DILC map. We use the newly- 
developed representation 1DFR of the spherical harmonic coefficients for cross and peak 
correlation. The 1DFR has an advantage in that the m = modes do not alter the 
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Figure 7. The top panel shows the cross-correlation between the foreground residuals 
Ri and the foregrounds F\ and the other two are the foreground contamination 
parameter q in units of cr^ mb . Middle panel is the contamination ratio Q/cr£ mb against 
the cross correlation coefficients X^ for DILC and the foregrounds (from Equation 
(5.11)). The curves are for X C J = —0.8 (the top) downwards with increment 0.2 to 
Xf = 0.8. The thick line is for zero correlation: X C J — between CMB and the 

foregrounds. In the bottom panel, we take the mean value of Xf^ from top panel of 
Figure 4 as input and perform 10 6 Monte Carlo simulations for the CMB and the X^ . 
The yellow, light blue and dark blue areas illustrate the contamination ratio e^/cr^ mb 
at 1, 2 and 3 a thresholds, respectively. 
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morphology, and when using the standard definition of cross correlation of the ID curves, 
they are intrinsically excluded in the calculation. 

We have tested the DILC map for 2 < £ < 10, although any diagnostics of 
correlation can be applied to all £ modes for both the WMAP and the upcoming 
Planck data. We find that eight out of nine modes are negatively-correlated, for which 
the probability is only 0.0176 if they are really independent. We also find the local 
extrema of the 1DFR curves are correlated, particularly for £ = 2 and 6 modes, with 
significance level less than 0.04. Both analyses have indicated that the DILC quadrupole 
is contaminated with foregrounds (with cross correlation coefficient X < —0.8 and 
extremum statistics at significance level a < 0.04). We then used the cross correlation 
coefficients to explore the level of the foreground contamination. 

The axiom that the CMB has no correlation with the foregrounds provides the most 
fundamental criterion that should be applied as a foreground cleaning check before any 
Gaussianity tests are performed and cosmologists get carried away by suggestions of new 
physics. We suggest that for the upcoming Planck mission, one should not ignore the 
kind of information we have presented here, before any scientific conclusion is reached. 

One final remark for our results is that the DILC map is derived from the concept 
of minimum variance optimization under the a priori assumption that CMB and the 
foregrounds are statistically independent. That we have only one universe, hence there 
is inevitably accidental correlation shall introduce some error in the estimation of the 
weighting coefficients, and consequently introduce residuals in the resultant map. Even if 
the correlation between the CMB and the foregrounds happens to be zero, the variance 
of the foreground residuals is minimum, but not necessarily zero. Although we can 
extract the angular power spectrum as close as possible to that of the CMB using such 
methods, the signal obtained in general still contains residuals. 
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